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We present exact models for an antiferromagnetic S = 1 spin chain describing trimerization as 
well as for an antiferromagnetic S — 3/2 spin chain describing tetramerization. These models can 
be seen as generalizations of the Majumdar-Ghosh model. For both models, we provide a local 
Hamiltonian and its exact three- or four-fold degenerate ground state wavef unctions, respectively. 
We numerically confirm the validity of both models using exact diagonalization and discuss the low 
lying excitations. 
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I. INTRODUCTION 

Spin 1 antiferromagnetic chains have been the subject 
of extensive study in the 1980s and 90s. Many studies 
were motivated by Haldane's identification of the 0(3) 
non-linear sigma model as the effective low-energy the- 
ory for spin chains |l| . He pointed out that antiferromag- 
netic spin chains with integer spin representations pos- 
sess a finite energy gap in the excitation spectrum and 
that the ground state correlations exhibit an exponential 
decay. Haldane's conjecture was substantiated through a 
rigorous theorem by Affleck and Lieb [2( and through the 
AKLT model dUll on theoretical site, and through the 
observation of the Haldane gap in spin 1 chains on exper- 
imental site [a, LD, la, lSl, llfll • Haldane's prediction has also 
been confirmed in detail by numerical studies [Til. IT3. [l3| . 

Recently, spin 1 chains have seen a considerable re- 
newal of interest. Experimental realizations using po- 
lar molecules stored in optical lattices have been pro- 
posed with the spin represented by a single- valence 
electron of a heteronuclear molecule. Moreover, it ap- 
pears possible that the one-parameter family of Hamil- 
tonians 



N 



He = J2 cosO SiS i+1 + sm9(S t S l+1 ) 2 
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can be engineered in optical lattices using cold spin 
1 bosonic particles with antiferromagnetic interactions, 
such as 23 Na, for arbitrary values of 9 [TEI. Hfl|. 

The phase diagram of the model ([JJ as a function of 9 
has been investigated by numerous authors ( e. q. fTjL [ill 
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references therein) and is by now well understood. 

The point 9 = on the circle shown in Fig. [1] the 
antiferromagnetic Heisenberg point, is embedded in the 
so-called Haldane phase (— 7r/4 < 9 < 7r/4) which is 
characterized by a unique ground state, exponentially de- 
caying correlations, and a gap between the ground state 
and the excited states. The Haldane phase includes at 
#vbs = arctan(l/3) the valence bond solid (VBS) or 
AKLT model. The AKLT Hamiltonian shares the most 
properties of the isotropic Heisenberg Hamiltonian but, 



in contrast to the isotropic Heisenberg model, possesses 
a ground state which can be written out explicitly. 

Above the Haldane phase in Fig. fJJ, there is a crit- 
ical phase (7r/4 < 9 < tt/2) with spin nematic corre- 
lations [3l|. The phase transition at $uls = 7r/4 was 
proposed to be of Kosterlitz-Thouless type 0, [3(| • At 
the transition point, the Hamiltonian UJ reduces to the 
Uimin-Lai-Sutherland (ULS) model [HS [H[ which 
exhibits explicit SU(3) symmetry. The ULS model is a 
sum of permutation operators and exactly solvable via 
the nested Bethe ansatz. 

At 9 = tt/2, the Hamiltonian fTJ becomes ferromag- 
netic with gapless excitations. It reaches the ferromag- 
netic Heisenberg point at 9 = ±7r and undergoes a first 
order phase transition to a dimerized phase at 9 = — 37r/4 
where JTJ is again SU(3) symmetric and has a highly de- 
generate ground state [36[ . Close to this point there was 
a long-standing discussion regarding the possible exis- 
tence of a small spin nematic phase. Recently, this was 
ruled out by numerical and analytical arguments (3ll . |3^ | . 
In the dimerized phase (— 37r/4 < 9 < — tt/4), the exci- 
tations are gaped. At the Takhtajan-Babudjan point 
#tb = — 7r/4, the gap closes and the model is again ex- 
actly solvable via the nested Bethe ansatz [38|, [3f| , has 
gapless excitations, and a unique ground state. Finally, 
the phase transition to the Haldane phase at 9 = —n /A 
is of second order [l7], E(| ■ 

Fath and Solyom [22|, H3 observed in 1991 a period 
tripling in the spectrum of |JJ in the critical phase 
(7r/4 < 9 < 7r/2). The observation of three soft modes in 
their numerical studies caused a controversy whether or 
not it is a trimer phase [U, [42], EH . Subsequent numerical 
studies concluded that there is no trimer phase [44]]. Re- 
cently, it was found that the dominant correlations in this 
phase are not of singlet, but of spin nematic (quadrupo- 
lar) character [3l| . 

Regardless of the spin nematic phase, S = 1 models 
yielding trimerization can be constructed with a spin in- 
teraction beyond the nearest-neighbor case. Solyom and 
Zittartz [45| presented such a model with four-site inter- 
action. In this model, the trimer singlets are nested and 
each trimer-singlet is placed on three non-neighboring 
sites, say on sites 3i — 2, 3i, and 3i + 2 (i = 1, . . . , N/3). 
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FIG. 1: Phase diagram of Hamiltonian |T} as a function of I 



Most recently, Corboz et at [4jg| investigated numerically 
the bilinear-biquadratic Heisenberg model with nearest- 
and next nearest neighbor interactions: 
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U = Ji cos e S,S l+1 + sin 6 (SiS i+1 ) 



i=l 
N 
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-J a ^ cos0SiS i+ 2 + sin 8 (S l S l+2 f 



i=l 



For certain values of the ratio J%/ J\, they found in a 
small region around (9tjls a trimerized phase in which 
the ground state becomes three-fold degenerate as they 
approach the thermodynamic limit. This can be seen 
as the analogy of the dimer phase in the S = 1/2 J±- 
J 2 -model for J 2 /Ji > 0.2411 [H H E|. Recently, 
trimerization and tetramerization was also discussed by 
Lecheminant and Totsuka [EH within a field theoret- 
ical approach where they considered a self-dual SU(n) 
Sine-Gordon model. 

In the first part of this article, we discuss an exact 
model for trimerization in an S — 1 antiferromagnet. 
The model is exact for finite (and hence trivially also 
for infinite) chains. The ground state consists of simple 
trimer singlets (each of which is placed on three con- 
secutive sites) and is three-fold degenerate. Since the 
three ground states can be written as valence bond solids, 
the model is in the spirit of the Majumdar-Ghosh (MG) 
model [52|, [53| for S = 1/2, where dimers occupy pairs 
of neighboring sites. We present a local Hamiltonian for 
which the trimer states are the exact zero-energy ground 
states, and identify numerically the elementary triplet 
excitation. It is a soliton consisting of two antisymmetri- 
cally coupled spins on adjacent sites. In the second part, 
we discuss the possible generalization to higher spin and 
introduce an exact model for the tetramerized (or, alter- 
natively, quadrumerized) phase in a spin 5 = 3/2 anti- 
ferromagnet. We numerically confirm the validity of the 
model using exact diagonalization. We further discuss 



low-lying excitations and give an outlook for tetramer- 
ization in other models. 



II. MODEL FOR TRIMERIZATION 

We consider a chain with N = 3/i sites (fi integer) and 
impose periodic boundary conditions (PBCs). To write 
the trimer ground states as products of trimer singlets, 
we define an operator T[i,j,k] which creates a trimer 
singlet on sites i, j, and k as follows: 

T[i,j,k] |0) = -^(|+,0,-) + |0 ) - +) + |-,+,0) 



|+,-,0)-|-,0,+)-|0, +,-)), 



(3) 



where | + ,0, — ) = cj + cj cj. _ | ) etc., with the usual 
fermionic creation operators cj , a = +,0,—. Using 
equation Q the three trimer ground states are given by 



= Y[T[3i+l,3i + 2,3i + 3}\0) , 
n-x 

\ip 2 ) = n T[3z + 2, 3i + 3, 3i + 4] |0) , (4) 

i=0 

l^s) = ]J T[3i + 3,3i + 4,3i + 5] |0) . 



Note that the trimer ground states ((3]) break transla- 
tional symmetry spontaneously, while they are invariant 
under translations by three lattice spacings. To illustrate 
the states, we show one of them in Fig. H]by connecting 
the sites belonging to the trimer singlets by arrows. 

To find an exact parent Hamiltonian for the trimer 
states, interactions including first-, second-, and third- 
nearest neighbors are required. More generally speaking, 
we conjecture that for the construction of an exact Hamil- 
tonian for "n-merization" consisting of local projection 
operators, interactions involving n+1 neighboring sites 
are required. Assuming interactions involving less than 
n + 1 sites, additional non-n-merized states might occur 
which will be annihilated by the Hamiltonian. We further 
stress that SU(3n) symmetry is a crucial requirement for 
trimerization (even though we are unable to prove this 
statement rigorously). 
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FIG. 2: Illustration of one of the trimer ground states on a 
chain with N — 3fi sites, where three neighboring spins are 
antisymmetrically coupled. 
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For convenience we introduce the auxiliary operator 
. It acts on four consecutive sites, 

X^ = £ S 3 S r + (S 3 S r ) 2 , (5) 

j;j'=i 
3<j' 

where Sj is the usual spin 1 operator. In terms of these 
the Hamiltonian we propose is given by 

N 

tt = £(Af>-6)(*W-4). (6) 

i=l 

As in the case of the MG Hamiltonian, equation ^ is a 
sum over projection operators. The trimer ground states 
will be annihilated by these operators, 

(Af>-6)(Af>-4)hM=0 

for v = 1, 2, or 3, and thus TL = 0. 

The key to the trimer phase is the explicit SU(3) sym- 
metry at the ULS point in Fig.[T] The symmetry emerges 
because the bilinear-biquadratic 5 = 1 Heisenberg in- 
teraction becomes proportional to the SU(3) symmetric 
permutation operator if the bilinear and the biquadratic 
term appear with the same coefficients (as it happens 
at the ULS point). The SU(3) symmetric permutation 

operator V^l fulfils 

^ ) |+,0) = |0,+), v®\+-) = \- +), 
V[f |0, -) = I-, 0) , V[f \+, +) = | + , +) etc., 

(3) 

i.e., V a l permutes the spins on site a and (5. The per- 
mutation operator and the spin 1 operators are related 
by 

Vf & = S a S p + (S a S p ) 2 - 1. (7) 

On the other hand, the permutation operator is related 
to the SU(3) spin operators J a by 

Vfp = 2J a J fj + 1 -. (8) 

The SU(3) generators at each lattice site a are defined 
as 

J « = 2 ^ C aaKa> C aa>> a = 1, . . . , 8, 

CT,O-'=+,0,- 

where A a are the SU(3) Gell-Mann matrices (see e.g. 
(EH). Equations j7|) and allows us to define a spin 
5 = 1 model which is simultaneously an SU(3) model. 
The ground states of the model are given by the trimer 
products (j4|). Note that §5§ is up to an additive and 
multiplicative constant equal to the Casimir of the total 
SU(3) spin on four consecutive sites; the Hamiltonian © 



corresponds to equation (12) in Ref. [54|. We now ex- 
plain the explicit construction for the SU(3) Hamiltonian 
which is equivalent to (|6|) and which describes trimeriza- 
tion. (A detailed discussion of this Hamiltonian for the 
SU(3) system can be found in [H|].) 

We consider a spin chain with a fundamental represen- 
tation (1,0) of SU(3) on each lattice site. As mentioned 
above, for exact trimerization interactions involving four 
neighboring sites are required. To find the relevant SU(3) 
representations appearing on four consecutive sites, we 
couple four fundamental representations, (1, 0) (1, 0) 
(1,0)0(1,0) = 3- (1,0)02- (0,2)03- (2,1)0 (4,0). 
In the ground state, only the representations (1,0) and 
(0, 2) are present. This can be seen by considering four 
neighboring sites in the ground state configuration (see 
Fig. [2]): either three spins are antisymmetrically coupled 
toasinglet {i.e., A[ (1, 0)® (1, 0)® (1, 0) ]<g> (1, 0) = (1,0)) 
or pairs of spins are antisymmetrically coupled to anti- 
fundamental representations (0, 1) (i.e., A[(l, 0)0(1, O)]0 
A[(l, 0) (1, 0)] = (0, 1) (0, 1) = (1, 0) © (0, 2)). Hence, 
the projection operator onto the subspace (2, 1) and (4, 0) 
applied to one of the trimer ground states must be zero, 

Prrm (i,j,k,l) |^) = Prjn (*,i,M) Vl>u) = (9) 

for v = 1,2, or 3 if i, j, k, and / label four consecutive 
sites. The Young table aux 1 1 1 1 1 corresponds to the rep- 
resentation (4,0) and [j- 1 - 1 to the representation (2,1). 
The conditions © singles out the states (|1]) uniquely 
as ground states, which enables us to write the parent 
Hamiltonian as 

W= E (Prrm M) + Prm (10) 
(ijkl) \ ~ ) 

The brackets (•) indicate summation over four neighbor- 
ing sites along the chain. We can rewrite this Hamilto- 
nian in terms of SU(3) spin operators, such that it annihi- 
lates the states which carry exclusively the representation 
(1,0) or (0,2) on four consecutive sites: 




where we have used that the eigenvalues of the squared 
total spin on four neighboring sites, ^4 S it es j is f m case of 
the representation (1,0) and in case of the representa- 
tion (0,2). The eigenvalues of the Casimir operators for 
the representations (2,1) or (4,0) are both larger than 
4p The Hamiltonian (fTTj) annihilates the ground states 
while the other states end up with a positive energy. 
Finally, we replace the SU(3) spin operators via ([5]) and 
d?]) by spin 5=1 operators to obtain ([6]). 

We have confirmed our predictions with exact diag- 
onalization (ED) for chains with N = 9 and N = 12 
sites. In both cases we found precisely three zero-energy 
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FIG. 3: Spectrum of the spin 1 model ((6]) for trimerization 
on a chain with N = 18 sites. The zero energy ground states 
(labeled by thick lines in the spectrum) are at k — 0, 4p 
in the Brioullin zone. 



ground states. These ground states are in the Brioullin 
zone located at k — 0, reflecting the fact that 

the system is translationally invariant under translations 
by three lattice spacings. We have plotted the spectrum 
E(k) in Fig. [3] for N = 18 sites where we used Lanczos 
routine for diagonalization. 

Even though we cannot write down any of the excited 
states exactly, we are able to elaborate on key properties 
like the quantum numbers involved. To create an exci- 
tation, we inevitably have to break a trimer, i.e., cre- 
ate a domain wall between the degenerate ground states. 
There are, however, two different types of domain walls 
(see Fig.0](a) and (b)) which correspond to different ex- 
citations. 

As explained above, our model exhibits an SU(3) sym- 
metry, and we may view the trimer singlets JU) as SU(3) 
singlets. Breaking such an SU(3) singlet yields either an 
individual SU(3) spin with fundamental representation 
(1,0) or two antisymmetrically coupled spins on adja- 
cent sites with the resulting anti-fundamental represen- 
tation (0,1), i.e., -4[(1,0) <g> (1,0)] = (0,1). Note that 
both (1,0) and (0,1) are three-dimensional representa- 
tions. Since the only three-dimensional representation 
of SU(2) is the triplet, the excitation has to be a triplet 
regardless of which type it corresponds to. 

The first type consists of an individual spin S = 1, 
playing the role of a domain wall between two differ- 
ent trimer ground states, e.g. between and l^) (see 
Fig.0|(a)). 

The second type consists of two spins S — 1 on ad- 
jacent sites. They are coupled such that the inter- 
nal Hilbert space spanned by this excitation is three 
dimensional, i.e., \e + ) = ^= (|+, 0) - |0, +)), |e ) = 

->-!-+». and |e_) = ^ (|0, -) - |- 0)). 
This second type of triplet excitation is also as a domain 
wall between two different ground states, e.g. between 
and ^3) (see Fig.@](b)). Since either type of do- 
main wall could in principle decay into two domain walls 



of the other type, only one type can be an eigenstate of 
the Hamiltonian ^ above. 

In Ref. (54j we have this question numerically investi- 
gated for the corresponding SU(3) model. We have found 
evidence that the domain wall placed on adjacent sites 
as shown in Fig.|5](b) is the elementary excitation. As 
pointed out above, this domain wall corresponds to an ex- 
citation transforming according to the anti-fundamental 
representation (0, 1) of the SU(3) model. This is what 
we expect as the spinon excitations for antiferromagnetic 
SU(3) spin chains generally transform according to the 
representation (0,1). For the Haldane-Shastry model, 
this was shown by exp licit construction of the exact one 
spinon eigenstates [55|,[5i|. On more general grounds, the 
low-energy behavior of an SU(3) spin chain with fun- 
damental representation is described by the SU(3)/- = i 
Wess-Zumino-Novikov-Witten (WZNW) model. The 
elementary excitations in this model transform likewise 
accor ding to the representation (0, 1) under SU(3) rota- 
tions daiEi]. 

Note that the triplet excitations are gaped, as it costs 
a finite energy to break a trimer. In Fig. [5] we have plot- 
ted the gap size for the finite chains versus the inverse 
system length. Fitting data points for N = 6, 9, 12, 15, 
and 18 yields a huge gap A = 4.56 ± 0.04 for the ther- 
modynamic limit. This gap, however, is not the Haldane 
gap exhibited by the antiferromagnetic S — 1 Heisenberg 
chain or the AKLT model, which we understand as due 
to a confinement force between spinons. This is consis- 
tent, as the SU(3) Heisenberg model does not display a 
Haldane gap [2l]54|. The static spin-spin correlations of 
our model © decay abruptly, as adjacent trimer singlets 
are uncorrelated. 

So far, we have used the mappings and ([5]) to con- 
struct a model of an antiferromagnetic S — 1 spin chain 
with an SU(3) symmetry. As a possible generalization, 
we might look for models of spin S antiferromagnets with 
an SU(2S , + 1) symmetry, i.e., models in which n = 2S+1 
spins might form an n-mer rather than three neighboring 
spins a trimer. Corboz et al. [46[ advocated that such an 
n-merized phase should appear in an appropriate J\—J% 
model provided the ratio J\j J% exceeds a certain critical 
value. Note that this is not in contradiction to our pre- 
vious statement that n + 1-site interactions are required 
for tetramerization. Corboz et al. proposed an n-merized 
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FIG. 4: (a) Possible triplet excitation as a domain wall on 
a chain with 3/i + 1 sites, (b) Possible triplet excitation as 
a domain wall on a chain with 3/j, + 2 sites. The excitations 
consists of two antisymmetrically coupled spins on adjacent 
sites. 
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FIG. 5: The gap size of a chain with certain length N is 
plotted vs. the inverse chain length. Data points for N = 6, 
9, 12, 15, and 18 yield the fitted value of A = 4.56 ± 0.04. 



In analogy to the MG and the trimer model ([6]), a 
Hamiltonian which annihilates the tetramer states can 
be written as a sum over projection operators. For con- 
venience we introduce an auxiliary operator acting 
on five neighboring sites 



i+4 

E" 

j<i' 



,(4) 
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phase which becomes exact in the thermodynamic limit 
only. Our statement is restricted to models describing 
exact n-merization. 

The generalization of our construction is based on the 
mapping 



? (n) 

a/3 



E^ l) 

p=0 



(12) 



where the S a are spin S operators and "P^a with n = 
25 + 1 is the SU(n) symmetric permutation operator, ex- 
ists in general. A method to determine the constants al" 
has been developed by Kennedy [5{|; Itoi and Kato [301 ] 
obtained them explicitly up to S = 2. We now use these 
results to introduce a spin 3/2 parent Hamiltonian which 
describes tetramcrization. 



III. MODEL FOR TETRAMERIZATION 

We consider a spin 5 = 3/2 antiferromagnetic chain 
with N = AfjL sites (/i integer) and PBCs. The operator 
Q[i,j,k,l] creates a tetramer singlet on sites i,j,k, and 



Q[i,j,k,l]\0) = 
1 



E J - - + 1 

Q,/3, 7 ,5=7r(|, I, -!,-§) 



C i,a C i+l,l3 C i+2,j C i+3,8 I ® i 



(13) 



The sum in (fT3"|) extends over all 24 permutations tt of 
the four states ||>. = c[, |0), ||>. = c[, |0), |-|>. = 

C t_ 4 |0), and |-|>. = 4_ f |( 
ground states are given by 



The four tetramer 



^> = II 0[4»+v,4i+l+i/,4i+2+i',4»+3+i/] |0) (14) 



i=0 



with ;y = 1,2,3, or 4. One of these ground states is 
illustrated in Fig. [S] 



Note that V^l is the SU(4) symmetric permutation op- 
erator whereas the Sj are spin 3/2 operators. In terms 
of the auxiliary operator the parent Hamiltonian is given 
by 



^quaclr 



E(^ 5) 



5)(i; 



(5) 



(16) 



This Hamiltonian is positive semi-definite, exact by con- 
struction, and annihilates the states 02 



X. 



(5) 



(->) 



= 0, 



which implies W quadr |^> = for v = 1,2,3, or 4. It 
represents an exact model for tetramerization in a spin 
3/2 antiferromagnet. 

We have confirmed our predictions with exact diago- 
nalization (ED) for chains with N = 8 and N — 12 sites. 
The spectrum of the 12 site chain is shown in Fig. [71 
where the four zero-energy ground states are located at 
k = 0,5, tt, and 4^, reflecting that the ground states are 
invariant under translations by four lattice spacings. 

Since the Hamiltonian (|16[) is similar to ^ , we con- 
tend ourselves with a brief discussion. Interactions be- 
tween five neighboring sites are required to ensure anni- 
hilation of each tetramer singlet and, hence, annihilation 
of the states (fT4")l . In order to find the correct projection 
operators, we couple five SU(4) spins with fundamental 
representation: (1, 0, 0) ® (1, 0, 0) <g> (1, 0, 0) ® (1, 0, 0) ® 
(1,0,0) = 4- (1,0, 0)85- (0,1, 1)06- (2,0, 1)©5- (1,2,0)© 
4 • (3, 1, 0) © (5, 0, 0). In the tetramer ground states, only 
the representations (1, 0, 0) and (0, 1, 1) are present. The 
projection operators onto the subspaces in which the rep- 
resentations (1,0,0) and (0,1,1) are absent hence anni- 
hilate the states (fT4")l . Since these states are the only 
states which contain representations (1, 0, 0) and (0, 1, 1) 



FIG. 6: Illustration of a tetramer state on a chain with JV = 
4/i sites, where four neighboring spins are antisymmetrically 
coupled. 
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FIG. 7: Spectrum of the spin 3/2 model (|16[) for tetramer- 
ization on a chain with iV = 12 sites. The zero energy 
ground states (labeled by thick lines in the spectrum) are 
at k — 0, 5, 7r, 4^ in the Brioullin zone. 



on five consecutive sites, a Hamiltonian in terms of pro- 
jection operators may be written as 

W quadr = E (^FFP {i,j,k,l,m) + Prrn (U.M.m) 

{ijklm) \ U 



+P \ (i,j,k,l,m) + P rjrn (i,j,k,l,m) 



where (•) indicates summation over five neighboring sites 
along the chain. We then replace the projection opera- 
tors by SU(4) spin operators and rewrite the SU(4) spin 
operators by spin S = 3/2 operators using the mapping 
(fT2| . This yields the Hamiltonian (fT6|) . 

Whereas the model for trimerization exhibited two 
possible candidates for the low-lying excitation, the 
model for tetramerization exhibit three candidates. 
When breaking a tetramer singlet which we might view as 
an SU(4) singlet, one obtains either an individual SU(4) 
spin with fundamental representation (1, 0, 0), or two an- 
tisymmetrically coupled spins on adjacent sites with rep- 
resentation (0,1,0), or three antisymmetrically coupled 
spins on consecutive sites with anti-fundamental repre- 
sentation (0,0,1). Whereas the representations (1,0,0) 
and (0, 0, 1) are four-dimensional, the representation 
(0, 1,0) is six-dimensional. The excitations with funda- 
mental and anti-fundamental representation correspond 
in case of the spin S = 3/2 chain to the quadruplet ex- 
citation carrying spin S = 3/2. However, if one couples 
two spin S = 3/2, the decomposition into irreducible 
representations does not contain a six-dimensional rep- 
resentation: !®| = 0ffil®2©3. We may hence discard 
the representation (0, 1,0) from the list of candidates. 

In order to motivate that the quadruplet excitation 
placed on three consecutive sites is the elementary ex- 
citation, we consider an SU(4) spin chain. For the 
SU(n) Haldane-Shastry model (HSM) (and particularly 



for the SU(4) HSM), the lowest lying excitation car- 
ries the anti-fundamental representation [56j . Moreover, 
the low-energy behavior of SU(4) spin chains with fun- 
damental representation is described by the SU(4)fc = i 
WZNW model where the lowest lying excitation carries 
the quantum number (0,0,1), i.e., it transforms under 
the anti-fundamental representation under SU(4) rota- 
tions [57], Hl| . This representation can be realized by an- 
tisymmetrizing three SU(4) spins with fundamental rep- 
resentation (1,0,0). It corresponds in the spin S=3/2 
model to the quadruplet excitation placed on three con- 
secutive sites (see Fig. [5]). This quadruplet excitation is 
gaped, as it costs a finite energy to break a tetramer 
singlet. 



FIG. 8: Quadruplet excitation serving as a domain wall be- 
tween two of the tetramer ground states on a spin S = § 
chain with 4fi + 3 sites. 



Most recently, the SU (4) J\-Ji model with fundamen- 
tal representation was investigated both within a level 
spectroscopy analysis of the ED data [6(| and within 
the density matrix renormalization group (DMRG) anal- 
ysis [sH. As proposed by Corboz et al. 4(|, the exis- 
tence of a tetramerized phase was verified for the regime 
Ji « J 2 @,[6l|. 



IV. CONCLUSION 

In conclusion, we propose exact models for trimeriza- 
tion and tetramerization in antiferromagnetic spin 1 and 
spin 3/2 chains, respectively. They can be seen as gen- 
eralizations of the Majumdar-Ghosh model. The mod- 
els consist of a local Hamiltonian involving four site or 
five site interactions, respectively, with a three-fold or 
four-fold degenerate ground state. The ground states are 
products of local trimer or tetramer singlets where each 
trimcr or tetramer is placed on three or four consecutive 
sites. We have numerically investigated the excitation 
spectrum and verified the validity of both models using 
exact diagonalization on finite chains. 

After this work was completed we became aware of 
another interesting generalization of the Majumdar- 
Ghosh model where three-anyon interactions in the con- 
text of a chain of Fibonacci anyons are investigated. 
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